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Abstract 

Using data collected with the BABAR detector between 1999 and 2002, we report the observation 
of 823±57 and 969±65 decays to D^*^D^*^K, where D^*^ and are fully reconstructed 
and where K is either a or a Kg decaying to vr^vr^. All 22 possible B decays to D^*^D'^*^K 
are reconstructed exclusively and the corresponding branching fractions or limits are determined. 
The preliminary branching fractions of the B^ and of the B^ to D^*^ D^*^ K are found to be 

13{B° D^*^D^*^K) = (4.3 ± 0.3(stot) ± 0.6(syst)) x lO^^ 
B{B+ D^*^D^*^K) = (3.5 ± 0.3(stat) ± Q.b{syst)) x lO-^. 
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1 Introduction 



The inconsistency between the measured b — > ccs rate and the rate of semileptonic B decays 
has been a long-standing problem in B physics. Until 1994, it was believed that the b — > ccs 
transition was dominated by decays B — > DgX, with some smaller contributions from decays to 
charmonium states and to charmed strange baryons. Therefore, the branching fraction b — > ccs was 
computed from the inclusive B ^ Dg X , B ^ (cc) X and B ^ ScX branching fractions, leading 
to B{b ccs) = (15.8 lb 2.8)% |jT|. Theoretical calculations are unable to simultaneously describe 
this low branching fraction and the semileptonic branching fraction of the B meson . 

As a possible explanation of this problem, it has been conjectured Q that B{b ccs) is in fact 
larger and that decays of the type B d''*^ D^*^ K {X) (where D^*^ can be either a D'^, D*^, 
or )0 could contribute significantly to the decay rate. This might also include possible decays 
to orbitally-excited Dg mesons, B — > D^*^ D**, followed by D** Z)^*) K. Experimental evidence 
in support of this picture has been published in the past few years. This evidence includes the 
measured branching fraction for wrong-sign D production, averaged over charged and neutral B 
mesons, by CLEO [S{B D X) = (7.9 it 2.2)%], and the observation of a small number of fully 

reconstructed decays B D^*^D'-*^K, both by CLEO g and ALEPH |]. More recently, BABAR 
and Belle ||8| have released some preliminary conference results on the evidence for transitions 
B^ D^*^^ D*~^ K~ with much larger data sets. 

B D^*^D'^*^K decays can occur through two different amplitudes: external W-emission 
amplitudes and internal W-emission amplitudes (also called color-suppressed amplitudes). Some 
decays proceed purely through one of these amplitudes while others can proceed through both. 
Fig. Q shows the possible types for charged and neutral B decays. In BABAR, the large data sets 
now available allow comprehensive investigations of these transitions. In the analysis described 
in this note, we present measurements of or limits on the branching fractions for all the possible 
B D^*^ D^*^ and B D^*^ D^*^ decay modes, using events in which both D mesons are 
fully reconstructed. 

2 The BaBAR detector and dataset 

The study reported here uses 75.9 fb"^ of data collected at the T{AS) resonance with the BABAR 
detector, corresponding to (82.3 it 0.9) x 10^ BB pairs. 

The BABAR detector is a large-acceptance solenoidal spectrometer (1.5 T) described in detail 
elsewhere |^]. The analysis described below makes use of charged track and vr*^ reconstruction 
and charged particle identification. Charged particle trajectories are measured by a 5-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH), which also provide ionisation 
measurements (dE/dx) used for particle identification. Photons and electrons are measured in the 
electromagnetic calorimeter (EMC), made of 6580 thallium-doped Csl crystals constructed in a 
non-projective barrel and forward endcap geometry. Charged K/tt separation up to 4GeV/c in 
momentum is provided by a detector of internally reflected Cherenkov light (DIRC), consisting of 
12 sectors of quartz bars that carry the Cherenkov light to an expansion volume filled with water 
and equipped with 10751 photomultiplier tubes. 



^Charge-conjugate reactions are implied throughout this note. 
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3 B candidate selection 



The and mesons are reconstructed in a sample of multihadron events for all the possible 
DDK modes, namely 5° ^ D^*^~ D^*^+ K'^ , D^*^^D^*)^K^ and B+ D^*^^ D^*)+ , 

D^*^^ D^*'^'^ , D^*^~ D^*'^'^ . mesons are reconstructed only from the decays Kg — > 7r'''7r~. 
To eliminate the background from continuum e~^e~ — > qq events, we require that the ratio of the 
second to zeroth Fox-Wolfram moments [10| be less that 0.45. 

The Kg candidates are reconstructed from two oppositely charged tracks consistent with coming 
from a common vertex and having an invariant mass within ±9 MeV/ of the nominal Kg mass. For 
most of the channels involving a Kg, we require that the Kg vertex is displaced from the interaction 
point for the event by at least 0.2 cm in the plane transverse to the beam axis direction. The vr'^ 
candidates are reconstructed from pairs of photons, each with an energy greater than 30 MeV, which 
are required to have a mass 115 < My-y < 150 MeV/ c^. The vr*^ from D*^ must have a momentum 
70 < p*(77) < 450MeV/c in the T(4S') frame, while the vr^ from D^ i^^vr+vr" must have an 
energy ^(vr") > 200 MeV. 

D* candidates are reconstructed in the decay modes D*'^ D^ir^, D*^ D^tt^ and D*^ 
Z)^7. A ±3(7 interval around the nominal mass difference AM = M{D*) - M{D°) is used to 
select D* mesons, where a is the measured mass resolution. For decays B^ — > D*~ D*~^ Kg and 
B+ D*-D*+K+, one of the D*^ is also allowed to decay to D^ir^. 

The D^ and D~^ mesons are reconstructed in the decay modes D^ — > K~'ir~^, K^tt^tt^, K^tt^-k^-k'^ 
and D~^ K~'k~^tt~^ , by selecting track combinations with invariant mass within ib2cr of the av- 
erage measured D mass. The average D mass and the D mass resolution a used in this selection 
are fitted from the data itself, using an inclusive sample of D decays. For modes involving two D^ 
mesons, at least one of them is required to decay to K^tt^ , except for the decay modes D*^ D*^ K^ , 
D*~ D*~^ K~^ and D*~ D^ K~^ , which have lower background. The K and n tracks are required to be 
well reconstructed in the tracking detectors and to originate from a common vertex. Charged kaon 
identification, based on the measured Cherenkov angle in the DIRC and the dE/dx measurements 
in the drift chamber and the vertex tracker, is used for most D decay modes, as well as for the K~^ 
from the B meson decay. 

B candidates are reconstructed from one D , one L>(*) and one K candidate. A mass constraint 
is applied to all the intermediate particles {D*^, D*^ , D^, D~^ , Kg, vr*^). Since the B mesons are 
produced via e~^e~ T(4S') BB, the energy of the B in the T{4S) frame is given by the beam 
energy E^^^^, which is known much more precisely than the energy of the B candidate. Therefore, 
to isolate the B meson signal, we use two kinematic variables: AE, the difference between the 
reconstructed energy of the B candidate and the beam energy in the center of mass frame, and 
mES) the beam energy substituted mass, defined as 



mES = \jE;^am-PB^ (l) 

where is the momentum of the reconstructed B in the T{4S) frame. Signal events have AE 
close to and tues close to the B meson mass, 5.729 GeV/c^. When several candidates are selected 
in an event, only the candidate with the lowest \AE\ value is considered ("best candidate"). Prom 
Monte Carlo studies, this algorithm is found to give the best reconstruction efficiency and the lower 
cross-feed between the different D^*^ D^*^K modes; it is found to introduce no bias on the signal 
extraction, since the latter is performed from the mss spectra only. However, in the Fig. to avoid 
the bias on AE inherent to the method, all the AE spectra are shown without this requirement. 
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4 Evidence for a signal in the sum of all B submodes 



The niEs and AE spectra of the selected events are shown in Fig. ^ for the sum of all the decay 
modes, separately for and B^. The AE spectra are shown for events in the signal region defined 
by 5.27 < mES < 5.29GeV/c^. Signal events appear in the peak near OMeV when reconstructed 
correctly, while the peak around — 160MeV is due to D DK or D D*K decays reconstructed as 
DDK or D DK, respectively. The ttt-es spectra for the signal region are shown for events with A^^ 
within ib2.5cTA£; of the central AE value for the signal. The resolution ctab is determined from 
the data and is equal to 9.9 MeV for events involving no D*^ and 11.3 MeV for events involving one 
D*^. For events with two D*^, the resolution is estimated from the Monte Carlo simulation to be 
13.8 MeV. A shift AE'ghift = (—5 it l)MeV of the AE central value for the signal is observed in 
the data. This shift is due to imperfect modeling of the charged K energy losses in the detector 
material and is accounted for in the analysis. As explained above, only the candidate with the 
lowest \AE — Ai^shiftl appears in the tties spectra in case of multiple candidates. Both the tties 
spectra for the AE signal region and the AE spectra show clear evidence of a signal. On the 
contrary, the rriES spectra for the background control region AE > 50 MeV show no evidence of 
any excess of events in the B signal region. For the ttt-es spectra, the combinatorial background is 



empirically described by the ARGUS function |ll|] 



dN 
druES 



f{mEs;A,C) = Ax mES 



m 



ES 



^beam 



X exp 



-CI 



m 



ES 



^beam , 



(2) 



where A is a normalisation factor. The function depends on a free parameter Q that is determined 
from a fit to the rriES spectra of the background control region. The number of combinatorial 
background events in the signal region is then estimated by normalizing the ARGUS function 
to the region 5.22 < rriEs < 5.27GeV/c^ in the projection containing the signal (Fig. §c,d) and 
extrapolating it to the signal region 5.27 < rriEs < 5.29GeV/c^. The fitted ARGUS functions are 
overlaid on the mES spectra of Fig. |2|. The average number of background events expected in the 
signal region is 1889ib24 for neutral B mesons and 2512ib27 for charged B mesons, while 2712 
and 3482 events are observed, giving an excess of 823±57 B^ and 969±65 B+ events in the signal 
region. 



5 Measurement of exclusive branching fractions 

In the following, the subscript k will be used to identify the different B D^*^D^*^K decay 
modes (i.e., 1^ D^K^, D*^ D^K^ , ...). The subscript i will be used to identify the different decay 
submodes of the DD pair (i.e., i =K'k x Ktt, Kit x Kittt^, Kit x KSir, ...). The subscript ik will 
therefore refer to the B mode k decaying into the DD submode i. 

The mES spectra obtained after a ib2.5fTA_E selection on (AE — AE'shift) for all the different 
D'-*^Di*)K modes are shown in Fig. ^ {B'^ decay modes) and Fig. ^ {B^ decay modes). The 
corresponding event yields, computed as explained below, are given in Table |l|. In Figs. ^ and ^ 
and in Table ||, for a given B decay mode the signals from the different DD decay submodes have 
been summed. However, to take advantage of the different signal-to-background ratio of the various 
submodes, the information from each submode is entered separately in a likelihood function used 
to calculate the B D^*^ D^*^ K branching fractions. As a first step, the ARGUS shape parameter 
of each submode, Qk, is determined from a fit to the mEs spectra of the background control region 



10 



/S.E > 50MeV. An ARGUS function with the shape parameter C fixed to this value is then fitted 
to the rriES distribution for the signal region \AE — A£'shift| < 2.5crA_B, excluding from the fit 
events with 5.27 < uies < 5.29GeV/c^. A value for the background normalization parameter Aik 
is calculated and the number of background events N^]^^ in the signal region for this submode is 
calculated as 

^ffc'= / fix;Ak,Qk)dx. (3) 

J5.27 

If submodes are used for a given mode, the branching fraction for that mode is then extracted 
by maximizing the following likelihood: 

i=l ^^^f^- 

where Nik is the observed number of events in the signal region and /Xj^ is the predicted number of 
events in the signal region. ^Uj^ is the sum of three contributions: 

• the predicted signal N^j^^ which is related to the (unknown) branching fraction Bk of decay 
mode A;, the reconstruction efficiency (ej^), the intermediate branching fractions and the 
number of BB events (N^-^) 

Nfk = BkX N^^ X e,k X Sf^; (5) 

• the number of combinatorial background events N^k^, determined as described above (Eq.^); 

• the peaking background Nfk'^^ from other B — > D^*'' D^*^K decay modes, calculated as 

= T.^ixNbbX ^'(^^ - X Bf^, (6) 

where e'{il ik) is the cross- feed matrix from B mode / to B mode k for the DD decay 
submodes i (the cross-feed between different DD decay submodes is found to be negligible). 
The only significant cross-feed is observed between decay modes where a fake D*^ replaces a 
true D*~^ or a true D^, for instance between D*~D^K~^ and D* D^K~^ , or between D* D^K~^ 
and D D*^K^. Therefore, these branching fractions are extracted with joint likelihood in 
Eq.|. 

The D* and D branching fractions used in the branching fraction calculation are summarized in 



Table ^ [12]. The selection efficiencies and the cross- feed matrices for each mode are obtained from 
a detailed Monte Carlo simulation, in which the detector response is modeled with the GEANT4 
program. B meson decays to DDK are generated with a three-body phase space model in the 
simulated event samples used for the efficiency calculation. For each decay submode, samples of 
about 15000 signal events have been produced. In addition, data are used whenever possible to 
determine detector performance: tracking efficiencies are determined by identifying tracks in the 
silicon vertex detector and measuring the fraction that is well reconstructed in the drift chamber; 
the kaon identification efficiency is estimated from a sample of D*'^ D^tt^, D^ K^tt^ decays; 
the 7 and vr*^ efficiencies are measured by comparing the ratio of events N{t~^ — > Vrh'^iT^)/N{T~^ — > 
) to the previously measured branching fractions |l^]. Typical efficiencies range from 
20%, for B° D^D'^K+ with both L>° mesons decaying to i^"7r+, to less than 1%, for B+ 
D*-D*+K+ {D*+ D^Tr+, D*- ifn- ) with mesons decaying to K-tt^-k^ or K--k+t:--k+. 
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Table 1: Number of events and branching fractions for each mode. The first error on the branching 
fraction is the statistical uncertainty and the second one is the systematic uncertainty. 







Total 






Branching 


90% C.L. 


B decay 


yield in 


Estimated 




fraction 


upper 


mode 




signal region 


background 


Excess 


(10-3) 


limit (10-3) 




B^ decays through external W- 


emission amplitudes 






D-D^K+ 


599 


479±12 


120±27 


1.7 ±0.3 ±0.3 






D-D*°K+ 


468 


337±10 


131±24 


4.6 ±0.7 ±0.7 




B^^ 


D*-D^K+ 


584 


399±11 


185±27 


3.lt°:^±0.4 




B^^ 


D*-D*^K+ 


289 


84±5 


205±18 


11.8 ± 1.0 ± 1.7 




B^ decays through external+internal W-emission amplitudes 


B^^ 


D-D+K^ 


26 


19±2 


7±5 


0.8+^:^ ±0.3 


< 1.7 


B^^ 


D*-D+K^+CC 


84 


34±3 


50±10 


6.5 ±1.2 ±1.0 




B^^ 


D*-D*+K^ 


116 


48±4 


68±11 


8.8^^:^ ± 1.3 




B 


decays through internal W-emission amplitudes 


B^^ 




175 


173±7 


2±15 


0.8 ±0.4 ±0.2 


< 1.4 






248 


225±8 


23±18 


1.7ti:3±0.7 


< 3.7 


B^^ 




123 


81±6 


42±13 


3.3i|i ± 1.4 


< 6.6 


S^ 


^ decays through external W- 


emission amplitudes 




B+ - 




367 


317±9 


50±21 


1.8 ±0.7 ±0.4 


< 2.8 


B+ - 




216 


175±7 


41±16 


4.lt}j±0.8 


< 6.1 


B+ - 


> D°D*+K^ 


77 


31±3 


46±9 


5.2^1:° ±0.7 




B+ - 




89 


43±4 


46±10 


7.8^2:? ± 1-4 




B^ decays through external+internal W-emission amplitudes 


B+ - 


> d'd^k+ 


627 


469±11 


158±27 


1.9 ±0.3 ±0.3 




B+ - 




552 


411±11 


141±26 


1.8l°i±0.4 


< 3.8 


B+ - 


> lfD*^K+ 


623 


402±11 


221±27 


4.7 ±0.7 ±0.7 




B+ - 




675 


468±15 


207±30 


5.3t\:l ± 1-2 






B- 


^ decays through internal W- 


emission amplitudes 




B+- 


> D-D+K+ 


64 


65±4 


-1±9 


0.0 ±0.3 ±0.1 


< 0.4 


B+ - 


. D^D*+K+ 


45 


39±4 


6±8 


0.2 ±0.2 ±0.1 


< 0.7 


B+ - 


^ D*-D+K+ 


64 


32±3 


32±9 


1.5 ±0.3 ±0.2 




B+ - 


> D*-D*+K+ 


83 


60±4 


23±10 


0.9 ±0.4 ±0.2 


< 1.8 



6 Systematic studies 

Due to the large number of and to the large multiplicities involved in the decays B 
D^*^D^*^K, the dominant systematic uncertainties come from our level of understanding of the 
charged kaon identification and of the charged particle tracking efficiencies. Both systematic uncer- 
tainties are estimated on a per track basis and are given in Table |3[ Other systematic uncertainties 
are due to uncertainties on the D and D* branching fractions, the ir^ reconstruction efficiencies, 
the D vertexing requirements, and the AE resolution used to define the signal box, as well as the 
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Table 2: Submode branching fractions used in the analysis |jT2|. The errors on the B{D0 — > 
K~-K^-K^) and B{D^ K~it^it^tt^) correlated to the error on B{D^ — > K~tt^) are indicated 
separately with the subscript Kt:. 



Mode 


B (%) 


K-TT+ 


3.80 ±0.09 




13.10 ± 0.84 ±0.31/^^ 




7.46 ± 0.30 ±0.18/^^ 


D+ K--K+TT+ 


9.1 ±0.6 


D*+ D^TT+ 


67.7 ±0.5 


D*+ D + TT^ 


30.7 ±0.5 




61.9 ±2.9 


D*o ^ 


38.1 ±2.9 


K^^ vr+vr- 


68.61 ±0.28 



uncertainty on the combinatorial background estimates, the statistical uncertainty on the efficiency 
due to the finite size of the Monte Carlo simulation samples and the uncertainty on the number 
of BB events in the data sample. The fractional systematic uncertainties on efficiencies and the 
branching fractions are summarized in Table ^. 

Possible decay model dependence of the efficiencies were also studied by generating decays 
^0 ^ D*-D+ and 5° ^ D*-D'+ {D+,D'+ D*^K+), where D^^ is the narrow {T = IMeV, 
M = 2.536 GeV/c^) orbitally excited 1+ state of the D** system and D'^ is a wide (T = 250 MeV, 
M = 2.560 GeV/c^) D** resonance. The efficiency for reconstructing these modes was compared to 
the efficiency found for decays B^ D*~ D*^ K'^ generated with a phase space model. We found 
no statistically significant difference in efficiencies; we assign a systematic uncertainty equal to the 
statistical error of 5%. 

7 Conclusions 

A preliminary measurement of the branching fractions for the 22 B ^ D^*^ D^*^K modes is given 
in Table [l|. For the channels for which S/\/B is smaller than 4, a 90% confidence level upper 
limit is also derived. (Here, B is the sum of the combinatorial background and of the cross-feed 
background from other D^*^ D^*^K modes and S = N — B, where N is the total yield in the signal 
region). This is the first time that a complete measurement of all the possible B D^*^ D^*'^K 
channels is performed. The measured branching fractions are in good agreement with earlier 
measurements made with smaller data sets for some of these modes ^, |^, For the decays 
proceeding through external W-emission or through the sum of external and internal W-emission 
amplitudes, the branching fractions B{B D* DK) and B{B — > DD*K) are found to be about 
twice the branching fraction B{B — > DDK). The branching fraction B{B D D*K) is found to 
be about 5 times larger than the branching fraction B{B — > DDK). No significant difference is 
observed between decays proceeding through external spectator amplitudes and decays proceeding 
through the sum of external and internal spectator amplitudes. 

After summing over all submodes, the preliminary branching fractions of the B^ and of the B~^ 



13 



Table 3: Fractional systematic uncertainties on efficiencies and branching fractions. 



Item 


Fractional uncertainty on efficiency or branching fraction 


Charged tracks reconstruction 


0.8% per track for good tracks in Drift Chamber 

1.2% per track for tracks without Drift Chamber requirements 


reconstruction 


2.5% per K^, added in quadrature to the track reconstruction error 


vr^ reconstruction 


5.1% per vr^^ 


T\ n 

7 from D*^ D^j 


5.1% per 7 (correlated with the vr^ systematic) 


identification 


2.5% per K± 


Vertex 

reconstruction 


1.3% per 2 track vertex 
3.1% per 3 track vertex 
5.7% per 4 track vertex 


a{AE) 


2% for modes with or 1 D*^ 
5% for modes with two D*^'s 


Monte Carlo statistics 


2% to 10% per DD submode (mode and submode dependent) 


Intermediate br. fraction 


see Table ^ 


Number of BB 


1.1% 


Decay model 


5% 



to D^*'d'^*^K are found to be 

B{B° D^*^D^*^K) = (4.3 ± 0.3(stat) ± 0.6(syst)) x 10"^ (7) 

B{B+ D^*^D^*^K) = (3.5 ± 0.3(stat) ± 0.5(syst)) x 10~^. (8) 

This study confirms that a significant fraction of the transitions h ccs proceeds through the 
decays B — > D^*^ D^*^ K . These decay modes account for about one half of the wrong-sign D 
production rate in B decays, B{B —^ D X) = (7.9 ± 2.2)% [Q]; however, because of the large 
statistical error on the latter measurement, it is not yet clear whether they saturate it. 

Future developments should include a search for D** resonant substructures in i? ^ d'^*'^ D^*^K 
decays, as well as a new high statistics measurement of the wrong sign D production in B decays 
and a search for decays B D^*'^D'^*^K* oi B ^ d'^*'^ D'^*^ K (mr) . 
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Figure 1: Top row: external W-emission amplitudes for the decays D * D^*^^K^ and 

D(*)-d(*^^K+. Second row: internal W-emission amplitudes for the decays B^ D^*^ D^*'>^K^ 

Bottom rows: external+internal W-emission amplitudes for the decays 
B+ ^ D^*^°d(*^^K+ and 5° ^ L>(*)-L>W+k0. 



16 



600 




0.1 0.2 

AE (GeV) 




0.1 0.2 

AE (GeV) 



600 




5.22 5.24 5.26 5.28 5.3 




(GeV/cO 



t^500 



IT) 

riOOO 

^500 




5.22 5.24 5.26 5.28 5.3 

(GeV/c) 



2000 
1500 
1000 
500 




- 1 1 1 1 1 1 1 


1 1 1 1 1 - 

(01 


" 1 1 1 1 1 1 


1 1 1 1 \ 1 " 



5.22 5.24 5.26 5.28 5.3 

nigs (GeX/c) 



Figure 2: The AE and mEs spectra (a,c,e) for the sum of all the S° ^ d'-*^D^*'^K modes and (b,d,f) 
for the sum of all the B+ D^*^D^*^K modes. (a,b): AE for mss > 5.27 GeV/ c^. (c,d): mES for 
AE < 2.5(T (signal box). (e,f): m-ES for AE > 50MeV (background control region). The curves 
superimposed to the m-Es spectra correspond to the ARGUS background fits described in the text 
and the shaded regions represent the background under the signal region 5.27 < rriEs < 5.29 GeV/c^. 
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Figure 3: The rriEs spectra of the ten — > Z) modes. For each mode, all the D decay 

submodes used in the analysis have been summed, except for plots where the D x D decay modes 
used appear explicitly. The curves correspond to the background fits described in the text and the 
shaded regions represent the background under the signal box. Upper row: pure external spectator 
^0 _^ £)(*)- j:)(*)Oj^+ decays. Middle row: external+internal decays 5° Bottom 

row: pure internal (color suppressed) decays D 
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Figure 4: The ttt-es spectra of the twelve D modes. For each mode, all the D decay 

submodes used in the analysis have been summed, except for plots where the D x D decay modes 
used appear explicitly. The curves correspond to the background fits described in the text and the 
shaded regions represent the background under the signal box. Upper row: pure external spectator 
decays B+ D^*^^ D^*^+ K'^ . Middle row: external+internal decays B+ d'-*^^ D^*^^ K+ . Bottom 
row: pure internal (color suppressed) decays D^*^~ D^*^~^ . 
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